Tuning Kondo physics in Graphene with gate voltage 
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We show theoretically that graphene, which exhibits a massless Dirac like spectrum for its elec- 
trons, can exhibit unconventional Kondo effect that can be tuned by an experimentally controllable 
applied gate voltage. We demonstrate the presence of a finite critical Kondo coupling strength in 
neutral graphene. We discuss the possibility of multichannel Kondo effect in this system which 
might lead to a non-Fermi liquid like ground state and provide a discussion of possible experimental 
realization of Kondo phenomenon in graphene. 

PACS numbers: 



Graphene, a two-dimensional single layer of graphite, 
has been recently fabricated by Novoselov et. al. In 
graphene, the energy bands touch the Fermi energy at 
six discrete points at the edges of the hexagonal Bril- 
louin zone. Out of these six Fermi points, only two are 
inequivalent; they are commonly referred to as K and K' 
points 2]. The quasiparticle excitations about these K 
and K' points obey linear Dirac-like energy dispersion. 
The presence of such Dirac-like quasiparticles is expected 
to lead to a number of unusual electronic properties in 
graphene including relativistic quantum hall effect with 
unusual structure of Hall plateaus Q. Recently, exper- 
imental observation of the unusual plateau structure of 
the Hall conductivity has confirmed this theoretical pre- 
diction [j. Further, the presence of Dirac-like quasipar- 
ticles in graphene provides us with an experimental test 
bed for Klein paradox , transmission resonance ^] and 
Lorenz-boost type phenomena. 

An extremely interesting phenomenon in conventional 
metal systems is the Kondo effect which occurs in the 
presence of dilute concentration of localized quantum 
spins coupled to the spin-degenerate Fermi sea of metal 
electrons [Tj . The impurity spin-electron interaction then 
results in perfect or partial screening of the impurity spin 
leading to an apparently divergent resistance, as one ap- 
proaches zero temperature. It also results in a sharp 
'Kondo Resonance' in electron spectral functions. Re- 
cent developments in quantum dots and nano devices 
have given new ways in which various theoretical results 
in Kondo physics, which are not easily testable otherwise, 
can be tested and confirmed experimentally Q. Most of 
the early studies in Kondo effect were carried on for con- 
ventional metallic systems with constant density of states 
(DOS) at the Fermi surface [Oj]. Some studies on Kondo 
effect in possible flux phases [lO[ , nodal quasiparticles in 
d-wave superconductors [ill . Luttinger liquids [l^ . and 
hexagonal Kondo lattice [13|, for which the DOS of the 
associated Fermions vanishes as some power law at the 
Fermi surface, has also been undertaken. However, al- 
though effect of non-magnetic impurities has been stud- 
ied lj|, there has been no theoretical study till date on 



the nature of Kondo effect in graphene. 

In this letter we present a large N analysis for a generic 
local moment coupled to Dirac electrons in graphene to 
show that Kondo effect in graphene is unconventional can 
be tuned by gate voltage. We demonstrate the presence 
of a finite critical Kondo coupling strength in neutral 
graphene. We point out that local moments in graphene 
can lead to non Fermi-liquid ground state via multi chan- 
nel Kondo effect. We also suggest possible experimental 
realization of such Kondo scatterers in graphene. 

The crucial requirement for occurrence of Kondo effect 
is that the embedded impurities should retain their mag- 
netic moment in the presence of conduction of electrons 
of graphene. We will not quantitatively address the prob- 
lem of local moment formation in the presence of Dirac 
sea of electrons in graphene in the present paper. We 
expect that large band width and small linearly vanish- 
ing density of states at the fermi level in graphene should 
make survival of impurity magnetic moment easier than 
in the conventional 3D metallic matrix. A qualitative 
estimate of the resultant Kondo coupling can be easily 
made considering hybridization of electrons in tt band 
in graphene with d orbitals of transition metals. Typi- 
cal hopping matrix elements for electrons in tt band is 
t ~ 2eV and effective Hubbard U in transition metals is 
8eV. So the Kondo exchange J ^ At^ /U , estimated via 
standard Schrieffer-Wolf transformation, can be as large 
as 2 eV which is close to one of the largest J ~ 2.5 eV 
for Mn in Zn. In the rest of this work, we shall therefore 
use the Kondo Hamiltonian [15] as our staring point. 

Our analysis begins with the Hamiltonian for non- 
interacting Dirac electron in graphene. In the presence 
of a gate voltage V , the Hamiltonian can be expressed in 
terms of electron annihilation operators ^^j-^-jQ 
lattice A{B) and Dirac point s — K,K with spin a =t, I 
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where sgn(s) = 1(— 1) for s = K{K'), vp is the Fermi 
velocity of graphene, and all repeated indices are summed 
over. In Eq. [T] and in rest of the work, we shall use 
an upper momentum cutoff = A/{hvp), where A ~ 
2eV corresponds to energy up to which the linear Dirac 
dispersion is valid, for all momenta integrals. 

Eq. [1] can be easily diagonalized to obtain the eigen- 
values and eigenfunctions of the Dirac electrons: E± — 
eV ± hvpk where k — (k^, ky) — (fc, 6) denote momenta 
in graphene and (m^"*", Hp"*^) — 1/a/2 (1, ± exp (isgn(s)0)). 
Following Ref. [10|], we now introduce the ^ fields, which 
represents low energy excitations with energies E±, and 
write 

na(k) = ^ui^cio = i/V2(e;a(k)+e„(k)), 

n^ik) = exp(*0)/y2(e; Jk) - e Jk)). (2) 

In what follows, we shall consider a single impurity 
to be centered around x = 0. Thus to obtain an ex- 
pression for the coupling term between the local mo- 
ment and the conduction electrons, we shall need to ob- 
tain an expression for ^(x = 0) = ^'(0). To this end, 
we expand the ^ fields in angular momentum channels 
5+a(k) = I]m=-oo whcre wc havc written 

k — {k,9). After some straightforward algebra, one ob- 
tains 

1 r'^" hrik 

n.io)-^l ^{^Uk)+ejk)). (3) 

Note that 'I'_b(0) receives contribution from m = ±1 
channel while for ^'^(0), the m = channel contributes. 
The Kondo coupling of the electrons with the impurity 
spin is given by 



Hk = 
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where g is the effective Kondo coupling for energy scales 
up to the cutoff A, S denotes the spin at the impurity 
site, r are the generators of the SU(iVc) spin group, and 
we have now generalized the fermions, in the spirit of 
large N analysis, to have Ns flavors (valley indices) Nf 
colors (sublattice indices) and Nc spin. For realistic sys- 
tems Nf — Nc ^ Ns = 2. Here we have chosen Kondo 
coupling g to be independent of sublattice and valley in- 
dices. This is not a necessary assumption. However, we 
shall avoid extension of our analysis to flavor and/or color 
dependent coupling term for simplicity. Also, the Dirac 
nature of the graphene conduction electrons necessitates 
the Kondo Hamiltonian to mix m = ± 1 and m = chan- 
nels (Eqs. [3] and HI). This is in complete contrast to the 
conventional Kondo systems where the Kondo coupling 
involves only m = angular momentum channel. 



The kinetic energy of the Dirac electrons can 
also be expressed in terms of the ^ fields: Hq = 

r E™=^oo E.,. {E+{k)iZ^^T^ + i?-(fc)C^^e^). 

Typically such a term involves all angular momenta 
channels. For our purpose here, it will be enough 
to consider the contribution from electrons in the 
m = 0,±1 channels which contribute to scattering from 
the impurity (Eqs. [3] and |4]). To make further analytical 
progress, we now unfold the range of momenta k from 
(0, oo) to (— cxD,oo) by defining the fields c^2)Q 

c'm.ik) = ^?t-^*^"(^»^(|fc|), fc>o, 

c'K2)a(fc) - fc<0, (5) 

so that one can express the ^ fields as '^A{B)a(^^ ~ 
I-oo ^V\kK{2)aik)- In terms of the c^^^-^^ fields, the ki- 
netic energy (in the m = 0,±1 channels) and the Kondo 
terms in the Hamiltonian can therefore be written as 



Ho ^ I ^ dk/{2^)Ekcl^cl 



/ V\k\VW\dkdk' 
-x{4l{k)T^^pctp{k')S'^), (6) 



where — eV -I- hvpk and summation over all repeated 
indices are assumed. 

Next we follow standard procedure fl6^ of representing 
the local spin by SU(A^c)Fermionic fields fa and write the 
partition function of the system in terms of the / and c 
fields 



Z 

So 
Si 



J VcVc^'VfVfVe e-'^'''\ S = So + Si + S2 
dk/{2^) (cfJ(fc,T)Go-icL(fc,T)), 
J/{4Tr^Nckl) dr / ^/\k\^/\k^\dkdk' 

Jo J-kc J-ka 

c!^{k,T)T2^ct^ik',T)f\iT)T^Js{r) 



0h ^.fce 

dr I 

J-k. 



where Gn 



dr [(/t (r) [hdr + e(r)] /„(r)) - e(r)g] , (7) 



hdr + Ek is the propagator for c fields. 



J — gNc/2 is the renormalized Kondo coupling, we 
have imposed the impurity site occupancy constraint 
Ea fafa — Q using a Lagrange multiplier field e(T). 

We now use the identity t^^t^s — ^c^asSp-y — SapS-ys 
[1^ and decouple using a Hubbard-Stratonovitch field 
(pi. In the large Nc limit one has S = So + S2 + S^ + S4, 
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where 

53 = dT 

Jo J-k 

54 = N,kl/J / dr</>r^(r)</>f(T) 

Jo 



(8) 



Note that at the saddle point level 



so that a non-zero value of 0^ indicates the Kondo phase. 
In what follows, we are going to look for the static saddle 
point solution with i/if (t) = and e(r) = cq 16]. In this 
case, it is easy to integrate out the c and / fields, and 
obtain an effective action in terms of (j)o and eo and one 
gets S' = S5 + Se with 

55 = -phN,TT[ln{th^r.- ^0 ~ NsNfroG'oituj,V)cPo)], 
Se = f3h (NsN.Nfkl |0o|' / J - eoQ) , (9) 

where Tr denotes Matsubara frequency sum as well as 
trace over all matrices and the Fermion Green function 
G'oiiPn, q) = G'q is given by [l3| 



G'r, 



-A 



2ll{hvF)' 



;{ipn-q)lii l/\ipn~q\ 



(10) 



where, in the last line we have switched to dimensionless 
variables p„ = hujn/A and q = eV/ A. 

One can now obtain the saddle point equations from 
Eq. El which are given by SS'/S<j)o = and dS'/Seo = 0. 
Using Eqs. [5] and [TOl one gets (after continuing to real 
frequencies and for T = 0) 



1/J 



Q/N, = l/(27r) / dpiy{p-iy-AoGo/2y 



(11) 



where we have defined the dimensionless variable Aq = 
NfN.l^oWiT^fi^vj,), p = hiu/A, Go = 2Tr{hvF)^G'o/A, 
I' = eo/A > 0, and have used the energy cutoff A for all 
frequency integrals. At the critical value of the coupling 
strength, putting v = and Aq = 0, we finally obtain 
the expression for Jc{q,T) 

Jciq,T) = J,(0) [l-2qln(l/g2)ln(fcBr/A)]"\l2) 

where the temperature ksT \s the infrared cutoff, Jc(0) = 
[nhvEkiy' / A = TT^A is the critical coupling in the ab- 
sence of the gate voltage, and we have omitted all sub- 
leading non-divergent term which are not important for 
our purpose. For = = g, we thus have, analogous 
to the Kondo effect in flux phase systems [10|, a finite 
critical Kondo coupling Jc(0) = tt^A ~ 20eV which is a 
consequence of vanishing density of states at the Fermi 
energy for Dirac electrons in graphene. Of course, the 
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FIG. 1: Sketch of the critical Kondo coupling Jc{q,T) as a 
function of temperature for several applied voltages q — eV/A. 
The Kondo phase exists for J > Jc- 



mean- field theory overestimates Jc- A quantitatively ac- 
curate estimate of Jc requires a more sophisticated anal- 
ysis which we have not attempted here. 

The presence of a gate voltage leads to a Fermi sur- 
face and consequently Jc{q,T) — > as T 0. For a 
given experimental coupling J < Jc(0) and temperature 
T, one can tunc the gate voltage to enter a Kondo phase. 
Fig. [1] which shows a plot of Jc{q,T) as a function of 
T for several gate voltages q illustrates this point. The 
temperature T*{q) below which the system enters the 
Kondo phase for a physical coupling J can be obtained 
using Jc{q,T*) — J which yields 

kgT* = Aexp[(l~ J,(0)/J)/(2(zln[l/(72])] (13) 

For a typical J ~ 2eV and voltage eV ~ 0.5eV, T* ~ 35K 
fv^ . We stress that even with overestimated Jc, physi- 
cally reasonable J leads to experimentally achievable T* 
for a wide range of experimentally tunable gate voltages. 

We now discuss the possible ground state in the Kondo 
phase qualitatively. In the absence of the gate voltage a 
finite Jc implies that the ground state will be non-Fermi 
liquid as also noted in Ref. [l3| for flux phase systems. 
In view of the large Jc estimated above, it might be hard 
to realize such a state in undoped graphene. However, 
in the presence of the gate voltage, if the impurity atom 
generates a spin half moment and the Kondo coupling is 
independent of the valley(flavor) index, we shall have a 
realization of two-channel Kondo effect in graphene ow- 
ing to the valley degeneracy of the Dirac electrons. This 
would again lead to overscreening and thus a non Fermi- 
liquid like ground state 0. The study of details of such 
a ground state necessitates an analysis beyond our large 
N mean-field theory. To our knowledge, such an analysis 
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has not been undertaken for Kondo systems with angu- 
lar momentum mixing. In this work, we shall be content 
with pointing out the possibility of such a multichannel 
Kondo effect in graphene and leave a more detailed anal- 
ysis as an open problem for future work. 

Next, we discuss experimental observability of the 
Kondo phenomena in graphene. The main problem in 
this respect is creation of local moment in graphene. 
There are several routes to solving this problem, i) Sub- 
stitution of a carbon atom by a transition metal atom. 
This might in principle frustrate the strong sp^ bond- 
ing and thus locally disturb the integrity of graphene 
atomic net. However, nature has found imaginative ways 
of incorporating transition metal atoms in p-7r bonded 
planar molecular systems such as porphyrin [isj . Sim- 
ilar transition metal atom incorporation in extended 
graphene, with the help of suitable bridging atoms, might 
be possible, ii) One can try chemisorption of transition 
metal atoms such as Fe on graphene surface through 
sp-d hybridization in a similar way as in intercalated 
graphite iii) It might be possible to chemically 

bond molecules or free radicals with magnetic moment 
on graphene surface as recently done with cobalt pthalo- 
cyanene (CoPc) molecule on AU(lll) surface [20|. This 
might result in a strong coupling between graphene and 
impurity atom leading to high Kondo temperatures as 
seen for CoPc on AU(lll) surface {Tk ^ 280K). iv) Re- 
cently ferromagnetic cobalt atom clusters with sub nano- 
meter size, deposited on carbon nanotube, have exhib- 
ited Kondo resonance 2l| . Similar clusters deposition in 
graphene might be a good candidate for realization of 
Kondo systems in graphene. v) From quantum chem- 
istry arguments, a carbon vacancy, or substitution of a 
carbon atom by a boron or nitrogen might lead to a spin- 
half local moment formation. In particular, it has been 
shown that generation of local defects by proton irradi- 
ation can create local moments in graphite 2^- Similar 
irradiation technique may also work for graphene. 

For spin one local moments and in the presence of 
sufficiently large voltage and low temperature, one can 
have a conventional Kondo effect in graphene. The 
Kondo temperature for this can be easily estimated using 
^bTk ~ -Dexp(— 1/pJ) where the band cutoff £) ~ lOeV, 
J ~ 2 — 3eV and DOS per site in graphene p ~ 1/20 per 
eV. This yield Tk ^ & - 150K. The estimated value of 
Tk has rather large variation due to exponential depen- 
dence on J . However, we note that Kondo effect due to 
Cobalt nano-particle in graphitic systems such as carbon 
nanotube leads to a high Tk ~ SOif which means that a 
large J may not be uncommon in these systems. 

Finally, we note that recent experiments have shown 
a striking conductance changes in carbon nanotubes and 
graphene, to the extent of being able to detect single 
paramagnetic spin-half NO2 molecule [23|. This has been 
ascribed to conductance increase arising from hole doping 
(one electron transfer from graphene to NO2). Although 



Kondo effect can also lead to conductance changes, in 
view of the fact that a similar effect has been also seen 
for diamagnetic NH3 molecules, the physics in these ex- 
periments is likely to be that of charge transfer and not 
local moment formation. 

In conclusion, we have shown that the Kondo effect 
in graphene is unconventional and can be tuned by an 
applied gate voltage. We have shown that it is possi- 
ble to have multichannel Kondo effect in graphene and 
discussed experimental possibilities of its realization. 

After submission of the first version of this manuscript 
in the arXiv ( |arXiv :0705.0257, vl), we became aware of 
Ref. 2J] with similar conclusion regarding existence of 
finite critical Kondo coupling in neutral graphene. 
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